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Abstract. In this talk, I will show how to determine the biasing factor 
b from the high-order moments of galaxies. The determination is based 
on the analytical modeling of primordial peaks and virialized halos and 
is independent of the currently unknown density parameter fio and other 
cosmological parameters. The observed high-oder moments of the APM 
galaxies require that the biasing factor b be very close to 1, i.e. the optical 
galaxies are an unbiased tracer of the underlying mass distribution (on 
quasilinear scale). The theoretical argument can be easily generalized to the 
three-point correlation function and the bispectrum both of which can used 
as further observational tests to the important conclusion of b ~ 1 drawn 
from the high-order moments. Finally I present our preliminary results of 
the three-point correlation functions for the Las Campanas Redshift Survey. 



1. Introduction 

A fundamental problem in cosmology is to understand how the spatial dis- 
tribution of galaxies (and of galaxy clusters) is related to that of the under- 
lying mass. The possibility for a galaxy bias (i.e. galaxies do not faithfully 
trace the underlying matter in spatial distribution) is advocated both by 
observations and by theories. It is therefore very important to determine 
the relation between the distribution of galaxies and that of dark matter 
(DM). The distribution of DM cannot be directly measured by telescopes 
because of its very nature. People attempted to understand the DM dis- 
tribution by studying the dynamics of galaxies since the peculiar motions 
of galaxies are closely related to the distribution of DM. Unfortunately, 
these studies can only yield a combination of O^os where Qo is the den- 
sity parameter of the Universe and as is the rms tophat density contrast 
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of radius 8h~ Mpc. Moreover, there are considerable discrepancies in the 
^o' 6 C8 value among different measurements. In this talk, I will show how 
to determine the biasing factor b from the high-order moments of galaxies, 
which is based on the analytical modeling of primordial peaks and virial- 
ized halos and is independent of the currently unknown density parameter 
and other cosmological parameters. The method is distinct from the 
related work of Gaztanaga & Frieman (1994) which is based on a simple 
non-physical parameterization of the bias relation. 

2. Measuring the biasing factor b 

The peak model and the halo model are generally used to prescribe the 
formation of galaxies. Recently, we have analytically derived the high-order 
moments for the primordial peaks and virialized halos (Mo, Jing, & White 
1997) based on the method of Mo &; White (1996). The high order moments 
of the peaks or the halos, e.g. the skewness Ss(R), the kurtosis S^R) and 
Ss(R), are functions of the biasing factor 6, the shape of the linear power 
spectrum, and the formation time Zf only. Our analytical model has been 
tested against a set of high-resolution cosmological N-body simulations, 
with a conclusion that the model works very well in the quasilinear regime. 

In Figure 1, I present our predictions for the high-order moments Sj t9 
(j = 3, 4, 5) of peaks and halos as a function of the linear bias parameter 
b. In this plot, a CDM-like power spectrum with T = 0.2 is assumed. This 
spectrum is known to be consistent with the angular correlation function of 
the APM galaxies (Efstathiou, Sutherland & Maddox 1990). The currently 
best estimates of the high-order moments for optical galaxies are those of 
Gatzanaga (1994) based on the APM survey, which are plotted on the figure 
with error-bars. Recent observations of galaxies at high redshift indicate 
that bulk of galaxies is formed at Zf > 1. Figures 1 then implies that the 
optical galaxies are little biased relative to the distribution of underlying 
mass (i.e. b s=s 1). It is interesting to note that the observations of S3, S4 
and S5 for the APM galaxies are all consistent with the analytical model if 
b « 1. 

A comparison with the observed moments of galaxy clusters is also 
interesting. Clusters of galaxies are generally believed to be halos formed 
recently (zf = 0) and have a bias factor b « 3 ~ 5. Our model predicts the 
skewness £3^ « 2 and the kurtosis S^h ~ 8 for such halos which are well 
consistent with the observations for the Abell clusters (Jing & Zhang 1989) 
and for the APM clusters (Gatzanaga, Croft & Dalton 1995). 

3. Further observational tests: the three-point correlation func- 
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Figure 1. Model predictions for the high order moments Sj t „ (j = 3, 4 and 5) of halos 
(left panels) and peaks (right panels) at a radius r — 10 h~ Mpc as a function of the 
linear bias parameter b. Each curve shows result for a given formation time Zf = 0, 0.5, 
1, 2 and 4 (curves from bottom up). The horizontal lines show the values of Sj for the 
mass density field calculated from quasilinear theory, whereas the data points (plotted 
arbitrarily at 6 = 2) show the observational results for APM galaxies (Gaztanaga 1994) 



tion and the bispectrum 

With a simple extension of the argument of Mo et al. (1997), we can derive 
the bispectrum and the three-point correlation function for the DM halos 
and the primordial peaks. In this extension we assume the bias relation by 
derived Mo et al. (1997) holds at all scales. We have tested this extension 
against a set of high-resolution N-body simulations with 1.7 x 10 particles, 
and have found that the extension is very successful (Jing 1997, in prepa- 
ration). Figure 2 illustrates such a test, where we compare the normalized 
bispectrum of N-body dark matter halos (symbols) with the predictions 
of the extension. The extension predicts that the optical galaxies should 
closely follow the underlying dark matter in the bispectrum and in the 
three-point correlation function, if the biasing factor of the galaxies is about 
1 as the high-order moments indicate (§2) and the bulk of the galaxies is 
formed before redshift 1. Both the amplitude and the shape dependence 
of the three-point correlation function and the bispectrum are potentially 
useful for constraining the bias factor and the linear power spectrum of the 
cosmological models (Fry 1994; Jing & Borner 1997; Matsubara & Suto 
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Figure 2. The normalized bispectrum Q(ki, 1(2, 8) of virliazed halos estimated from the 
scale-free simulations with the power spectrum index n and 256 3 particles (symbols), 
compared with the analytical predictions discussed in the text. The halos of mass M (in 
unit of the particle mass) are formed at the scale factor a\ and the clustering analysis is 
made at a 2 . 

1994; Matarrese et al. 1997). 

4. Preliminary results from the LCRS survey 

Here I present the preliminary results of our determination (Jing & Borner 
1997) of the three-point correlation function for the Las Campanas Redshift 
Survey. With this large sample we are able to determine the three-point 
correlation function of galaxies much more accurately than any previous 
works (e.g. Efstathiou & Jedrzejewski 1984). The three-point correlation 
function is not only a powerful test to the bias problem as described above 
but also a generally important statistic to describe spatial distributions and 
dynamics of galaxies (Peebles 1980; Suto 1993). 

We analysed two types of three-point correlations: the redshift space 
three-point correlation function ((s, u, v) and the projected three-point cor- 
relation function U(r p ,u,v). Our results are presented in Figure 3. The 
figure shows that the redshift three-point correlation function is approxi- 
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Figure 3. The normalized three-point correlation function of galaxies in the Las Cam- 
panas Redshift Survey. The left panel shows the result in the redshift space and the right 
panel shows the projected one. The r, u, v convention of Peebles (1980) is adopted in 
the plot, and different symbols are for different values of u. More detailed account of the 
figure will be given in Jing & Borner (1997) 



mately hierarchical with the amplitude Q(ied) w 0.6 (see e.g. Jing & Borner 
1997 for its definition). There is a tendency that Q increases with v espe- 
cially for large triangles with s > 5 /i _1 Mpc, qualitatively in agreement with 
the model predictions for the mass three-point correlation function (Jing & 
Borner 1997). Similar properties are found for the projected one. However 
the projected Q(proj) is about 1.5 for very small scales r p < 0.5/i _1 Mpc 
and decreases with r p for larger r p . If the three-point correlation function 
in real space is hierarchical, the projected one should also be hierarchical 
with Q(P r °j) = Q- Out results then support that the three-point corre- 
lation function on small scales is hierarchical with Q ~ 1.5, consistent 
with the conclusions of Peebles (1980) based on 2-D angular catalogues. 
The hierarchical form however seems not to hold for scales much larger 
than 0.5/i _1 Mpc because Q(proj) decreases with r p and increase with v 
for r p > 0.5/i _1 Mpc. A quantitative comparison between the observations 
and the CDM models is in progress, which will take account of the redshift 
distortion effect and the projection effect (Jing & Borner 1997). 
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5. Conclusions 

In this talk I have presented our analytical model for calculating the high- 
order moments of peaks and halos in the quasilinear regime. With this 
model, the observed skewness and kurtosis of the Abell and the APM clus- 
ters are interpreted as a consequence of these clusters being recently formed 
massive halos. The high-oder moments of the APM galaxies are consistent 
with the model predictions if the galaxies are an unbiased tracer of the 
underlying mass distribution (on quasilinear scale). This model can even 
be extended to the normal three-point correlation function and the bispec- 
trum (cf. Catelan et al. 1997), and thus our important conclusion that 6^1 
can be further tested with these functions. Finally I presented the prelim- 
inary results of our determination of the three-point correlation functions 
for galaxies in the Las Campanas Redshift Survey. 
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